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We present a theoretical and experimental investigation of homonuclear ionizing collisions of laser- 
cooled metastable (2 ^Si) helium atoms, considering both the fermionic ^He and bosonic *He isotope. 
The theoretical description combines quantum threshold behavior, Wigner's spin-conservation rule 
and quantum statistical symmetry requirements in a single-channel model, complementing a more 
complete close-coupling theory that has been reported for collisions of metastable *He atoms. The 
model is supported with measurements (in the absence of light fields) of ionization rates in magneto- 
optically trapped samples, that contain about 3 x 10* atoms of a single isotope. The ionization 
rates are determined from measurements of trap loss due to light-assisted collisions combined with 
comparative measurements of the ion production rate in the absence and presence of trapping light. 
Theory and experiment show good agreement. 



I. INTRODUCTION 

Soon after the first demonstration of laser cooling of 
neutral atoms, it was recognized that collisions between 
atoms have a profound effect on the physics of laser- 
cooled atomic gases Q. Light-assisted collisions and 
other inelastic collision processes constitute loss channels 
that limit the attainable atomic densities in laser-cooled 
samples. The magneto-optical trap (mot) was found to 
be a versatile tool for accurate investigation of these in- 
elastic collisions that occur at collision energies of 
10~^ eV in samples with temperatures around 1 mK. 
Measuring loss rates from the trap, the cross section for 
inelastic collisions can be determined with great preci- 
sion Experiments as well as theoretical work, almost 
exclusively carried out in alkali systems, have provided 
a wealth of information on cold collisions 0, that play 
a critical role in many research areas that have emerged 
since the advent of laser cooling Q . 

Laser-cooled helium atoms in the metastable 2 '^Si 
state, denoted He*, provide unique opportunities to 
study cold ionizing collisions. 
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where pi stands for Penning ionization and Ai for associa- 
tive ionization. The simple structure of the helium atom 
. ^ allows for theoretical exercises that are relatively uncom- 
, plicated, while experiments profit by the possibility of 
^ • direct detection of ions with charged-particle detectors. 
5^ \ Furthermore, large numbers of atoms (~ 3 x 10^), either 
■^He* or "^He*, can be confined in a mot so that differ- 
ences between collisions of fermionic ('^He*) and bosonic 
atoms ("^He*) can be investigated. 

In this article, a theoretical and experimental inves- 
tigation of homonuclear ionizing collisions of ^He* and 
''He* atoms is presented, with particular regard to iso- 
topic differences in ionizing collisions. The theoretical 
model is a single-channel calculation of the ionization 
cross section and rate coefficient for ionizing collisions, 
that can be applied to both isotopes. In the exper- 
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iments, magneto-optically trapped samples of ^He* or 
■^He* atoms are investigated and ionization rates are mea- 
sured with a micro-channel plate (mcp) detector. The 
measurements are performed in the absence of trapping 
light, as optical excitation by near-resonant light alters 
the dynamics of cold coUisi' 
enhancing the ionization rate 

Combining theory and experiment, the present work 
extends previous studies and may be used to unravel 
inconsistencies in reported experimental results Ell- 
The transparent theoretical model complements close- 
coupling calculations that have been performed for colli- 
sions of *He* atoms El El E3- Other studies of ^He* 
collisions provide little detail |1, or focus exclusively 
on collisions of spin-polarized atoms (to investigate the 
feasibility of Bose- Einstein condensation) El ^3 j while 
another theoretical study, that applies to both '^He* and 
''He* collisions EIIj is based on questionable assump- 
tions. The extensive description of the experimental pro- 
cedure and exhaustive explanation of the data analysis 
procedure provides clear insight into the origin of exper- 
imental ionization rate coefficients and can be used to 
evaluate other experimental studies of ionizing collisions 
of He* atoms H El El 113 • 

The article is organized as follows. Section ITTl presents 
the theoretical model for homonuclear ionizing collisions 
of laser-cooled helium isotopes. First a simplified expres- 
sion for the ionization cross section is derived (Sec. lII A|l . 
Using an effective molecular potential (Sec. llllJ)l . partial 
wave ionization cross sections are derived from numeri- 
cal solutions of the one-dimensional Schrodinger equation 
(Sec. Ill C|) and the corresponding partial wave ionization 
rate coefficients are calculated (Sec. Ill D|) . Quantum sta- 
tistical symmetry requirements are taken into account 
(Sec. Ill E|) to derive the rate coefficient for samples of 
either ^He* or ''He* atoms with a given population of 
magnetic sublevels. Finally, the results are compared 
with other theoretical results (Sec. IIIFl) . Section ^O] de- 
scribes the measurement of the ionization rate coefficients 
in laser-cooled samples of '^He* or ''He* atoms. Subdivid- 
ing loss mechanisms into ionizing and non-ionizing contri- 
butions, and distinguishing between linear and quadratic 
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trap loss, an overview of trap loss mechanisms occurring 
in our He* samples is presented and estimates of ion- 
ization and trap loss rates are derived fSec. It is 
shown that ionization rates can be deduced from trap loss 
measurements if the contribution of ionizing mechanisms 
to trap loss is determined. Measurements of trap loss 
due to light-assisted collisions are presented fSec. IIIIBll . 
Comparative ion production rate measurements in the 
absence and presence of trapping light are used to de- 
termine the ionization rate coefficients in the absence of 
light fields f Sec. 11110)1 . Finally, the results are compared 
with our theoretical predictions fSec. IIIID|l . 

Section IIVI presents a discussion of the theoretical 
and experimental results, as well as conclusions and 
prospects. 



II. THEORETICAL MODEL 

The ionizing collisions of Eqs. ^ and ^ are highly 
exothermic, as the internal energy of two He* atoms ex- 
ceeds the 24.6 eV ionization energy of the He atom by 
more than 15 eV. As differences between Pi and Ai are 
unimportant for the work presented here (the reaction 
mechanisms have been discussed in detail in Ref. H^l), 
we will not distinguish between them and use the term 
PI to denote both processes. 

The interaction that drives the auto-ionizing transi- 
tions of Eqs. and |(2J) is of an electrostatic nature 
|22| . so that it only induces transitions between molecu- 
lar states of equal total electronic spin. Therefore, ion- 
ization rates associated with the reactions of Eqs. Q and 
^ depend on the total spin states on the reactant and 
product side of the reaction formulas. For both reactions, 
the reactants carry an electronic spin of s = 1 and can 
form total spin states with S* = 0, 1 or 2, while the prod- 
ucts, carrying s — ^ (except for ground state helium, 
which carries no electronic spin), can only form states 
with S — or 1. Clearly, total electronic spin can only 
be conserved if 5 = or 1, and a pi reaction with 5 = 2 
would involve a violation of spin conservation. It has 
been shown 19, 20] that a very weak spin-dipole mag- 
netic interaction can induce spin flips and mediate Pi in 
collisions of He* atoms with S = 2. The corresponding 
ionization rate is four orders of magnitude smaller than 
those of collisions with S = and S* = 1, for which total 
electronic spin is conserved |23j . The strong suppression 
of PI by spin conservation is known as Wigner's spin- 
conservation rule 13] and has been observed for He* col- 
lisions in a gas discharge a laser-cooled sample 
and a Bose-Einstein condensate In heavier 

metastable rare-gas systems, Wigner's spin-conservation 
rule does not apply [29j, and polarized samples can be 
used for the investigation of quantum statistical effects 
in ionizing collisions [s^, |^ . 

At the mK temperatures of a laser-cooled sample of 
He* atoms, collisions occur at relative kinetic energies 
E = w 10"'' eV, where fi — m/2 is the re- 



duced mass of the colliding atoms (with m the mass of 
the He atom) and is the relative velocity between the 
atoms. As the de Broglie wavelength of atomic motion 
A = h/^v-[ ~ 250 ao (with h Planck's constant) is much 
larger than the typical scale of the interatomic potential, 
collisions are dominated by quantum threshold behavior 
HE nil ' In this case, the collision process can 
be described conveniently using the partial wave method 
|36j |: the ionization cross section, written as a sum of 
partial wave contributions, 

=^af"\ (3) 

is dominated by only a few partial waves £. For inelastic 
exothermic collisions, such as Penning ionizing collisions, 
the quantum threshold behavior of the £th partial cross 
section is given by ^ il^ <34 | 

af")cxA:2^-i if fc ^ 0. (4) 

Here, k ~ {2fiE/h'^y^^ is the wave vector of the asymp- 
totic relative motion of the colliding atoms, with h = 
h/2'K. In a sufficiently cold sample of He* atoms, the 
cross section for Penning ionizing collisions is dominated 
by the s-wave contribution ctq , which diverges as 1 /k if 
A: — *■ 32] . Elastic collisions have very different threshold 

properties: the cross section a^^'^'^^^ approaches a nonva- 
nishing constant, (j^^^**^) (x fc* and a'f}^^'' cx A:^, if A: — > 

For collisions of He* atoms, partial wave cross section 
cr^'°"'' can be derived from the solution of an effective 
one-dimensional potential scattering problem (Sees. 1X1"^ 
IIID|I . Restrictions imposed by the symmetry postulate 
on the partial waves that contribute to the cross section 
of Eq. Q can be taken into account thereafter (Sec. lIIE|l . 

A. Ionization cross section 

From a semi-classical point of view, two events can be 
distinguished in the process of a cold ionizing collision 
of two He* atoms: (1) elastic scattering of the atoms by 
the interaction potential V{R)^ with R the internuclear 
distance, and (2) Penning ionization that occurs when 
the two electron clouds start to overlap As colli- 

sion energies are small, the elastic scattering occurs at a 
relatively large internuclear distance R > 100 aq. For 
partial waves £ > 0, the radial wave function Uke{R) 
(cf. Eq. ©) is scattered by the centrifugal barriers, 
while scattering for £ = takes places at the inter- 
nuclear distance where the local de Broglie wavelength 
A(i?) = h/{2n[E - y(i?)]}i/2 becomes comparable to 
the size of the potential, i.e. dA{R)/dR « 1 [H H^]. 
As the electron clouds of both atoms start to overlap at 
small internuclear distance i? ~ 5 oq js^ , the elastic scat- 
tering process can be considered to precede the inelastic 
process of ionization. 
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FIG. 1: Ab initio potential energy curves in atomic units, 
calculated for "He* by Miiller et al. 



In the spirit of Ref. |3l| , we assume that the two sub- 
sequent processes can be treated separately. As pi is 
a strong inelastic exothermic process, we can write the 
ionization cross section for collisions with total electronic 
spin S as |3l| 



(2S+l)^('™) 



{21+1) 



(2S+l)p(t"") (2S+1) p('™) 



(5) 

where (2'S'+i)p(*"'^) jg ^j^g probability for the atoms to 
reach a small internuclear distance, and (2-S+i)p('°") jg 
the probability for ionization to occur at that place. 
As total spin S is strongly conserved during ionization, 
(2S+i)p('°") is very small (<C 1) for collisions that violate 
Wigner's spin-conservation rule. Here, we neglect ioniz- 
ing collisions with S = 2 and assume that 5p('°"' = 0. As 
ionization occurs with essentially unit probability for the 
other spin states (Miiller et al. j38| report an ionization 
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probability of 0.975), we set 

The calculation of cross sections (2S'+i)u'''""' jg reduced 
to the determination of partial wave tunneling probabil- 
ities (2S'-(-i)p('"")^ ^jjg energy dependence of the prob- 
gives rise to an energy dependent 
cr , that displays the quantum threshold behav- 



abilities (2S+i)pf"" 



ior of the inelastic collisions of Eq. To calculate 

(2S+i)p(*"")^ we need to consider the interaction poten- 
tial of the colliding atoms. 



obtained from Miiller et al. [331 . The curves are the re- 
sult of ab initio calculations in the Born-Oppenheimer 
approximation, where total electronic spin S* is a good 
quantum number. The possible values S — S — 1 
and S — 2 correspond to a singlet, triplet and quintet 
potential '^V{R) and ^V{R)), respectively. In 

Fig.n the curves are labeled in Hund's case (a) notation 
2S+i^±^^^ where A = \Ml\ with Ml the quantum num- 
ber of the projection of the total electronic orbital angu- 
lar momentum onto the internuclear axis of the molecule, 
g/u stands for gerade or ungerade, i.e. positive or neg- 
ative symmetry under inversion of all electronic coordi- 
nates of the molecule, and ± indicates positive or nega- 
tive symmetry under reflection through a plane includ- 
ing the internuclear axis . As both electrons of the 
He* atom are in s-states, colliding atom pairs can only 
have zero total orbital angular momentum, indicated by 
A = E. 

The potentials can be extended to large internu- 
clear distance using a calculation of the dispersion in- 
teraction of two He* atoms; for a multipole expan- 
sion —Cq/R^ — Cs/R^ — Cio/i?^°, dispersion coefficients 
Ce = 3276.680 a.u., Cg = 210566.55 a.u. and Cio = 
21786760 a.u. have been reported Ell. Here, we 
construct potentials vahd for R > 3.5 oq by fitting the 
short-range potential curves of Fig. ^ smoothly onto the 
long-range dispersion interaction around 20 ao by inter- 
polation using a cubic spline fitted to R^ x (^s+i) 

The elastic scattering for collisions with S = and 
5 = 1 is governed by potentials ^V{R) and '^V{R), re- 
spectively. Within the framework of the partial wave 
method, potential scattering by ^^^'^^^V{R) is described 
by the radial wave equation |42| | 



2^ di?2 



uki{R) 



2^i?2 



(2S+1) 



V{R) - E 



UkdR) = 0, (6) 



where i is the quantum number of the relative angu- 
lar momentum and Uke{r) is the radial wave function. 
Equation ^ can be interpreted as a one-dimensional 
Schrodinger equation (i? > 3.5 oq), describing the poten- 
tial scattering of a particle of mass fi by effective potential 



(2S+1) 



Vi{R) 



(23+1] 



V{R) 



2^i?2 ' 



(7) 



B. Effective potential 

As the helium atom, with only two electrons and a nu- 
cleus, has a relatively simple structure, interatomic po- 
tentials can be calculated with high accuracy. Figure 
shows the short-range part (3.5 ao < R < 14.0 ao) of the 
potential curves for two metastable 2 '^Si helium atoms 



where the second term is the well-known centrifugal po- 
tential. 

To calculate the probability (2S'+i)p^*"") of atom pairs 
to reach the distance where ionization occurs, we modify 
the effective potential curves to simulate the ionization 
process. We set the curves to a constant value for small 
internuclear distances, and extend the range of R to neg- 
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ative values 1311 



(25+1) 



R < Ro, 
R > i?o, 



(8) 



where Rq — 6.1 oq is chosen to be the location of the 
potential curve minimum. In this way, we avoid reflec- 
tions of the radial wave function from artificial features 
of the potential energy curve at Rq. Modeling the in- 
teratomic interaction by Ve{R), potential scattering is 
described by the one-dimensional Schrodinger equation 
{—CO < R < co) 



2^ di?2' 



Uki{R)+ i{R) - E um{R) = 0, (9) 



and ionizing collisions correspond to the transmission of 
the relative particle to i? < 0: for atoms that reach 
the region of small _R, where ionization takes place, the 
corresponding relative particle will propagate freely to 
R = —oo and never reflect back to i? > Rq. The dis- 
appearance of the particle to i? < results in a loss of 
probability flux. 

Figure El shows plots of V"o(^), ViiR), and V2{R) for 
homonuclear collisions of both "^He* and "^He* atoms, 
where we have used the atomic mass m = 3.01603 u for 
^He* and m = 4.0026 u for ''He*, with lu = 1822.89 a.u. 
For an atomic sample with a thermal velocity distribu- 
tion with temperature T, the mean collision energy is 
given by {E) = ^k^T; this relation is used to express the 
potentials in units of temperature. 

The barriers formed by centrifugal potentials with 
£ = 1 and 2 are five orders of magnitude smaller than 
the potential energy associated with the short-range at- 
traction of the colliding atoms. However, even the lowest 
barrier, with £ = 1, is large compared to the temperature 
of 1 mK that is typical of samples of laser-cooled He* 
atoms, and the barrier heights increase with increasing 

£. Therefore, probabihty (2S-Hi)p^*"") is small for £ = 1 
and decreases rapidly for larger values of £. 



C. Effective potential scattering problem 

The probability of transmission through a barrier in a 
potential can be calculated from the stationary states as- 
sociated with the potential . To calculate the proba- 
bility for transmission to i? < in potential ^'^^~^^'>Vi{R), 
we consider stationary states that correspond to the sum 
of an incident and a reflected wave for i? 3> 1, and a 
single transmitted wave for R < Rq- The transmission 
probability (2S+i)p(*"") ^q^j^ |-,g written as js^ 

(25+1) p(t-) ^ (10) 
in 

where Jjn and Jtr are the probability fluxes associated 
with incident and transmitted plane waves, respectively. 
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FIG. 2: Potentials '^^+^Vf(i?) with centrifugal barriers. The 
centrifugal interaction is smaller than the short-range attrac- 
tion by five orders of magnitude. 



Numerical methods are used to calculate (2S'+i)p^'""^^ f-Q^. 
£ — 0,1,2 and a range of collision energies E. 

For a typical collision energy E — 9.5 x 10^^ a.u. 
{E/^k-Q = 2.0 niK), the transmission probabilities are 



2 

(tun) 



0.66, 



(tu 



0.086 and ip^*"'') 



5.8x10- 



Clearly, reflection is almost complete in the case of d- 
wave scattering (as E is much smaller than the barrier 
height), while transmission is significant for p-wave scat- 
tering. In the case of s-wave scattering, there is con- 
siderable reflection, although a centrifugal barrier is ab- 
sent. Here, quantum reflection occurs due to the mis- 
match between the long asymptotic de Broglie wave and 
the rapidly oscillating wave in the region of small inter- 
nuclear separation. We checked that the resulting partial 
wave cross sections satisfy the quantum threshold behav- 
ior of Eq. (gl). 

To determine the dependence of the cross sections on 
the adapted short-range part of the potential, we have 
calculated the variation in the cross sections as a func- 
tion of Rq for various collision energies. If Rq is close 
to the location of the potential curve minimum at 6.1 oq. 
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the variations are smallest (less than 0.2%). Further- 
more, the difference between probabilities ip^'""^^ and 
at a given collision energy E is only a few per- 
cent, as potentials ^V{R) and ^V{R) differ very lit- 
tle in the region where elastic scattering takes place: 
\^V{R) - ^V{R)\/\^V{R)\ < 10-4 for R > 20ao. 

It has been shown in calculations that the ionization 
cross sections for ''He* are enhanced if the s-wave scat- 
tering length associated with the quintet potential is 
near a singularity ISj. The s-wave scattering length de- 
scribes elastic collisions in the low-temperature limit [s^ 
and shows a singularity (goes through ±oo) whenever a 
bound state is removed from the potential 0. From 
experiments, it has been determined that = 7.6 nm 
with an error of 0.6 nm |43l |. The scattering length is 
sufficiently far from the singularity to neglect enhance- 
ment of the ionization cross sections. For '^He* atoms, 
the s-wave scattering length for S* = 2 is predicted to be 
5.0 nm <^a < 6.0 nm ^14?|. This is also sufhciently small 
to neglect effects on the ionization cross sections. 



D. Ionization rate coefficient 

The ion production rate dNio^/dt in a magneto- 
optically trapped atomic sample of ■^He* or ^He* atoms 
can be expressed in terms of an ionization rate coefficient 
K (particle"^ cm'^/s). 



dNionjt) 

dt 



K n^{r,t)d\. 



(11) 



The rate coefficient depends on the temperature T of the 
sample and can be written as 

/■OC 

K{T)= / f7('°")(Wr) 4"^\wr) V.dVr, (12) 
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FIG. 3: Partial wave ionization rate coefficients ^ICt and un- 
polarized rate coefficient ii-f""?"') for ^He* and ''He* {^ICe 
coefficients are, within a percent, equal to corresponding ^K-e 
coefficients) . 



with 



TT (fcBr/A*)3/2 



Pr'M - ^/z 717^^ exp( , (13) 



kBT/fiJ 



the Maxwell-Boltzmann distribution for the relative ve- 
locity in the atomic sample under study, with {v^)^/'^ = 

{knT/^iY'^ = (2fcBr/m)i/2. 

Correspondingly, we can define partial wave ionization 
rate coefficients, 

poo 



(14) 

and the ionization rate coefficient associated with poten- 
tial (2S+i)v^(i?)^ 



(2S- 



-i)i^(T) = ^ 



(15) 



Figure |31 shows plots of the rate coefficients (calculated 
with numerical integration) for a temperature range from 



10 to 100 mK, for collisions of ■^Hc* or "^He* atoms. 
The energy of the centrifugal barriers for £ = 1 and £ = 2 
is indicated by dashed vertical lines. For temperatures 
T < 5 mK, the contribution of the d-wave becomes very 
small and can be ignored if an accuracy of 5% is sufficient. 

The quantum threshold behavior of the rate coeffi- 
cients is given by 



(2S+l)^ 



e cx 



(16) 



In particular, we find ^ constant if fc ^ 0, 

as the divergence of (Jq is canceled by v = hk/fx. The 
differences between the two isotopes are less than 50%. 

We have neglected the atomic hyperfine structure of 
the interatomic potentials for •^He*. As it is four orders 
of magnitude smaller than the attractive interaction at 
short range, its effect on (2S+i)^('°"' (2S+i)^^ jg neg- 
ligible. 
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E. Symmetrization of scattering states 

Although the interatomic interaction is almost iden- 
tical in the case of '^He* and ''He*, giving rise to par- 
tial wave contributions (2S+i)g-('°") g^j^j {2S+i)jq^ tYiaX are 
similar, the composition of the total ionization cross sec- 
tion or rate coefficient from these contributions is very 
different for the bosonic ('*He*) and fermionic isotope 
(•^He*). The symmetrization postulate requires that a 
scattering state describing a colliding pair of identical 
bosons has even symmetry under exchange of the atoms, 
while a state describing identical fermions has odd sym- 
metry [45I l46|. As a result, partial waves with improper 
symmetry do not contribute to the total cross section or 
rate coefficient, and are excluded from the summations 
in Eqs. Q and (fT^ . 



1. Bosonic He*: Symmetric states 

The total electronic spin S and relative angular mo- 
mentum I of two colliding '*He* atoms are, to good ap- 
proximation, conserved separately 19, 2Qi, and S and i 
can be considered good quantum numbers. Ignoring the 
radial part of the quantum states, a basis for atom pairs 
is given by 



[\{si)A{s2)B,SMsJ.me)'^ 



(17) 



For the moment, the atoms are assumed to be distin- 
guishable and are labeled with A and B. The atoms 
carry spins Si — 1 and S2 = 1, respectively, and the to- 
tal electronic spin S is the result of the addition of si 
and S2- Quantum numbers Ms and are associated 
with the projection onto the internuclear axis of the to- 
tal electronic spin and the relative angular momentum 
respectively. 

For a system of identical bosons, such as a pair of ^He* 
atoms, physical states are symmetric under exchange of 
the two atoms. Such states are obtained applying sym- 
metrizer S = (1 -I- Pi2)/\/2 to the basis vectors |3a.l47l|. 



\s1S2, SMsJmi) 

= ^[l + (-lf+'] \{s,)A{s2)B,SMs,lm,), (18) 



and normalizing the result if necessary. As the states 
differ from zero only ii S + £ \s even, we can conclude 
that it is not possible to construct states with the proper 
symmetry if = or 5* = 2 and £ is odd, or if = 1 
and £ is even. Consequently, the corresponding partial 
wave ionization rate coefficients are excluded from the 
summation of Eq. H15|) and we can write 



I even 



(19) 
(20) 



As total electronic spin is conserved, each long-range 
scattering state converts to a single short-range molecular 
state associated with a potential of Fig. Consequently, 
the strong suppression of ionization in the quintet poten- 
tial results in a rate coefficient ^^K 0. 



2. Fermionic He*: Antisymmetric states with hyperfine 
structure 

In a collision of two ^He* atoms, the total atomic an- 
gular momentum F = fi + £2 and the relative angular 
momentum of the two atoms £ are, to good approxima- 



tion, conserved separately. Here, f. 



+ h U - 1,2), 



the total angular momentum of an atom is the sum of 
its electronic spin and nuclear spin ij. 

Ignoring the radial part of the quantum states and as- 
suming for the moment that the colliding atoms are dis- 
tinguishable, a basis for atom pairs is given by 

{\isiHh)A{s2t2f2)B, FMpJme)} , (21) 

where atoms A and B have identical electronic spins. 



that add 



si = S2 = 1, and nuclear spins, ii = 12 — 2' 
up to /i for atom A and /2 for atom B. The latter two 
add up to the total atomic angular momentum F with 
projection onto the internuclear axis Mp. In magneto- 
optically trapped '^He* samples, all atoms occupy the 
lower / = I hyperfine level, so that /i = /2 = | and 
F can take values 0, 1, 2 and 3. Quantum numbers £ 
and mi are the angular momenta associated with the 
relative motion of the two atoms and its projection onto 
the internuclear axis, respectively. 

As the '^He atom is a fermion, the physical states de- 
scribing atom pairs are antisymmetric under exchange of 
the atoms. Applying antisymmetrizer A = (1 — Pi2)/V^ 
to the basis states ^3 , obtain 

\siiifi,S2i2f2,FMF,£mi) 

= ^ [1+ (-1)^-^] \{siiih)Ais2t2f2)B,FMF,£me). 

(22) 

These states, that must be normalized if necessary, dif- 
fer from zero only if _F — £ is even. Consequently, only 
even (odd) partial waves contribute to collisions with 
even (odd) F, and the total ionization rate coefficient 
can be written 



J2 w) if 0,2, 



odd 



K{F) = { (23) 
^ /C<.(F) ifi^ = l,3. 

e odd 

To express partial wave rate coefficient ICi^F) in terms 
of the rate coefficients ^^'^"'"^^/Cf that are associated with 
singlet and triplet potentials, we consider the interaction 
between two colliding '^He* atoms. 



7 



The cohision process of two laser-cooled ^He* atoms 
is controlled by the atomic hyperfine interaction and 
the various interatomic interactions [^S l4(l l4ll | that re- 
sult in potentials (^"^+^^^(7?). As a result of the hyperfine 
interaction, S is not a good quantum number for large 
internuclear distances, where atom pairs are character- 
ized by F. However, S* is a good quantum number for 
R < 30 ao, where the molecular interaction dominates 
and Wigner's spin-conservation rule applies. The evolu- 
tion of the quantum mechanical state from long to short 
internuclear distance is well approximated by a diabatic 
transition, as the absolute change in the coupling between 
quasi-molecular states, during one period of oscillation 
in the quasi-molecular system, is much larger than the 
absolute energy difference between the scattering states 
[49I Isof. Around R — 35 oq, the molecular interaction, 
increasing exponentially with decreasing R [Tsl l , becomes 
larger than the atomic hyperfine interaction, while the 
relative velocity has increased to 4 x 10~^ a.u. due to 
the attractive Van der Waals potential. Consequently, 
the asymptotic scattering state of a colliding atom pair 
remains unchanged when the atoms reach small internu- 
clear distance and the relation between rate coefficient 
ICe{F) and coefficients (^S+i)/^^ 

can be determined by 
expanding the corresponding scattering state onto the 
molecular states. 

The expansion coefficients only depend on the angular 
part of the quantum states involved. For the molecu- 
lar states, the angular part can be derived applying the 
antisymmetrizer to basis 

[\isiii)Ais2i2)B,SI,FMF,emi)^ (24) 

(and normalizing if necessary), where the atoms A and 
B are assumed distinguishable, / is the quantum number 
associated with the sum of the nuclear spins, I = ii -f 
12, and S is the quantum number associated with the 
sum of the electronic spins, S = Si -I- S2; the bases of 
Eqs. 1)21(1 and 1)24(1 are related through 9-j symbols [HTj . 
The resulting physical states 



\siii, 3212, SI, FMp, tmi 



V2 



[i + i-iy-'-'] 



X \{siii)A{s2i2)B,SI,FMF,eme) (25) 

are different from zero only if £ — S — I is even. The 
expansion of a scattering state onto the molecular states 
can be confined within the subspace defined by F and £, 
as F and £ can be considered good quantum numbers, 

8212/2, FMpJme) 

= J2asi{F) \siH,S2i2,SI,FMF,£me). (26) 

S.I 

Table presents the expansion coefficients asi{F) for 
scattering states with /i = /2 = |. 

The partial wave rate coefficient lCi{F) associated with 
scattering state \siiifi, Siiif2, FMp, img) can be writ- 



ten as a weighted sum of coefficients ^'^^'^^'^ICi asso- 
ciated with molecular states \siii, siii, SI , FMp, £1711) 
with weights \asi{F)\'^, 



/C,(F)=^|as/(i^)Px(25+i)/c,. (27) 

SJ 



It can be seen in Table that, in the case of F = 0, only 
singlet and triplet states are involved, so that 



/c,(o) = f (i/C,) -f i(3/C,). 



(28) 



For the F — 1 and F — 2 states, the diabatic transi- 
tion transforms the scattering states into a superposition 
of ionizing and non-ionizing molecular states. The con- 
tribution to the ionization rate coefficient from quintet 
states can be neglected, so that 



/c,(i) = i^(i/c,) K^M, 

-( 

3 * 



/C,(2) = f (i/C,) 



(29) 
(30) 



Finally, in the case of F = 3, partial rate coefficient ICi{3) 
is negligible, as only quintet states are involved. 



F. Ionization rate coefficient for trapped samples 

In a laser-cooled sample of He* atoms, collisions occur 
for all values of the total atomic angular momentum, F = 
0, 1, 2, 3 in case of ^He*, and 5 = 0, 1, 2 in case of "He*. 
The contribution of each collision channel depends on 
the distribution P,„ of magnetic substates in the sample, 
where m is the azimuthal quantum number of the atom, 
which can take on values irif = ~fi~5j^i| in case of 
■^He* and = —1,0,1 in case of ''He*. Using density 
operator js^ 



Pir) = X! Prn{r)Pnir)\m,n){m,n\ 



(31) 



m n<r 



to describe a statistical mixture of (properly sym- 
metrized) magnetic substate pairs |m,n), with m and 
n azimuthal quantum numbers, the ionization rate coef- 
ficient for the mixture (cf. Eq. I(ll() ') can be written as 

£ even 

+ J2 C^°dd ^ICi + %odd '/C,)] n(r) dV. (32) 

e odd 

where N is the number of trapped atoms, n{r) is the 
density distribution in the sample, and '•^"^"'"^•'feovcn/odd 
is the sum of the expectation values of the density op- 
erator for all ionizing molecular states with total spin 
S and even/odd parity. Explicit expressions for the co- 
efficients are given in Table ^ As the coefficients can 
be interpreted as projections of the statistical mixture 
onto subspaces of ionizing states, Eq. ((32(1 implies the 
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TABLE I; Expansion coefficients asi{F) — {s\i\, 8212, SI , FMf ,imt\s\i\f\, S2i2.f2, FMf ,imt) . Scattering state 

\siiifi,S2i2f2,FMF,£nie) is indicated by its values of F, wliile molecular state \sii-i, 8212, SI , FMf, Irrii) is denoted in the 
Hund's case (a) notation, ^'^+^E+^^ (7). 

iE+(/ = 0) iE+(/ = l) ^E+iI = 0) 3E+(/=1) 5j.+ (j^q) 5E+(/ = 1) 



F = 

F = 1 
F^2 

F = 3 



10 5 — 2l. 

27 V 9 V 27 



41 



TABLE II: Coefficients ''^'^''''^''6oven/odd from Eq. I|82^ . The coefficients are the sums of the expectation values of the density 
operator for all ionizing molecular states of given S and given parity. 

^He* *Hc* 

'feeven |(P-3/2J'3/2 + P-I/2P1/2) |(2P-lPl + Pq") 

'beven §(^-3/2^-1/2 + P- 3/2^1/2 + P-I/2P3/2 + P1/2P3/2) + 5 (P-3/2P3/2 + P-I/2P1/2) 

'&odd |(P-3/2Pl/2 + P-I/2P3/2) + ^(P-l/2P-l/2 + P1/2P1/2) + IP-3/2P3/2 + ^P-l/2Pl/2 

^bodd 3(P-3/2Pl/2 + P-I/2P3/2) + |(P-l/2P-l/2 + P1/2P1/2) + 5P-3/2P3/2 + l^P-l/2Pl/2 P-^lPo + P-lPl + PoPl 




FIG. 4: Ratio of unpolarized ionization rate coefficients 
^3Hc^°'V^4Ho''°'^ including partial wave contributions up to 



assumption that scattering states transform diabatically 
to a superposition of molecular states, when atoms move 
from large to small internuclear distance. 

The unpolarized ionization rate coefficients iir*^"'^P°'\ 
i.e. the rate coefficient for a laser-cooled sample of He* 
atoms where the magnetic substates of the atoms are 
uniformly populated, is obtained by setting P_i = Po = 
Pi = i for 4He* or P_3/2 = P-1/2 = P1/2 - P3/2 = \ 
for ^He*. For samples with a temperature around 1 mK, 
only s and p- waves have to be taken into account. For 
unpolarized '''He* samples, we obtain 



TABLE III: Calculated values of if (""p°i) and comparison be- 
tween various theoretical results. Theoretical results from 
Ref. 15] are extracted from Fig. 8 of that paper. 





Ref. 


T (mK) 


i4:(""P°''(cmVs) 


This work 


^He* 


LIS] 


0.001 


9.9 X 10~" 


7.7 X 10"" 




M 


0.5 


2.2 X 10~^° 


7.5 X 10"" 




[15j 


0.5 


8.6 X 10"" 


7.5 X 10"" 




m 


1 


7.3 X 10"" 


8.3 X 10"" 




fl5] 


1 


8.9 X 10"" 


8.3 X 10"" 






2 




9.9 X 10"" 


^He* 




0.5 

2 


1.0 X lO"'' 


2.0 X 10"^° 
1.8 X 10"^° 



where ^ICq and ^JCi are the partial wave ionization rate 
coefficients for ^He*, calculated in Sec. Ill Cl In the case of 
an unpolarized sample of '^He* atoms in the lower / = | 
hyperfine level, the rate coefficient is given by 

3He « 1^ m'iCo) + leiCo) + f (^/Ci) + f(3/Ci)l 



K 



(unpol) 
4Hc 



i [(i/Co) + 3(3/Ci)] , 



(33) 



K. 

(34) 

where ^ICq, '^JCq, ^/Ci and "^/Ci are the partial wave ion- 
ization rate coefficients for ^^He*. For both isotopes, the 
coefficients are shown in Fig.|31for temperatures between 
10 fiK and 100 mK. It can be seen in Fig.0] where the 
ratio of the unpolarized rate coefficients for the two iso- 
topes is displayed, that unpolarized ■^He* atoms have a 
larger rate coefficient than unpolarized "^He* atoms for 
temperatures between 10 fxK and 100 mK. 

In Table UTTl we compare the results of the theoretical 
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model presented here to results of other theoretical work 
IE ^1 • We see that our results agree well with the 
results of the detailed close-coupling theory of and 
the simpler calculation of but that there is a large 
discrepancy with the results of Kumakura and Morita 
This is not surprising, since their model does not 
account for quantum reflection for s-wave scattering. As 
we have shown in Sec. Ill CI quantum reflection is signifi- 
cant and we estimate that the omission of this effect leads 
to rate coefficients that are too large by factor of about 2. 
Moreover, Kumakura and Morita assume that the evo- 
lution of the scattering states during the collision of two 
■^He* atoms can be approximated by an adiabatic tran- 
sition and accordingly apply the noncrossing rule 52] to 
derive the number of ionization channels for each partial 
wave ionization rate coefficient. As shown in Sec. Ill El 
the system is well approximated by a diabatic transition 
and we estimate that the assumptions of an adiabatic 
transition leads to an unpolarized rate coefficient that is 
50% too large. This explains the difference between our 
results and those of Ref. 11]. As a final remark we note 
that the observed differences between our work and those 
of Refs. [3 and |0] may be due to theoretical uncertain- 
ties in the molecular potentials and in the form of the 
ionization widths used in the calculations. This leads to 
theoretical uncertainties of ~40%, as discussed in Ref. 

M. 



III. MEASUREMENT OF IONIZING 
COLLISIONS IN A MAGNETO-OPTICAL TRAP 

The experimental investigation of ionizing collisions is 
performed in a setup that can be used to confine large 
numbers (> 10*) of atoms in a mot. As the setup has 
been described previously ^ , only a brief outline is pre- 
sented here. The MOT is loaded from a coUimated and 
Zeeman slowed atomic beam, produced by a liquid nitro- 
gen cooled, DC-discharge source. The beam source can 
be operated with pure '^He* or pure "'He* gas. In the 
case of ^He*, the gas is recycled and purified using liquid 
nitrogen cooled molecular sieves. 

The laser light that is used for coUimation, Zeeman 
slowing and magneto-optical trapping has a wavelength 
of 1083 nm and is near resonant with the 2 '^Si(/ = §) ^ 
2 ^P2(/ = |) optical transition in the case of '^He*, and 
the 2 ^Si 2 transition in the case of ''He*. For 
both transitions, the natural linewidth r/27r = 1.62 MHz 
and the saturation intensity /gat = 0.167 mW/cm^ (cy- 
cling transition) . The light is generated by an ytterbium- 
doped fiber laser, that is frequency-stabilized to the 
laser cooling transition using saturated absorption spec- 
troscopy in an rf-discharge cell. Acousto-optic mod- 
ulators (aoms) are used to generate the slowing and 
trapping frequencies, which are detuned by —500 MHz 
and —40 MHz, respectively. The slowing beam is fo- 
cused onto the atomic beam source, while the trapping 
light is split up into six independent Gaussian beams 



with rms diameters of 27 mm and a total peak inten- 
sity 

-^peak 

59 mW/cm2 (Jpeak/^sat « 353). The 
magnetic quadrupole field of the MOT is generated by 
two anti-Helmholtz coils and has an axial field gradient 
dB/dz = 0.28 T/m. 

For the investigation of ionizing collisions, the trapped 
He* samples are studied with an absorption imaging sys- 
tem and two microchannel plate (mcp) detectors [50l] . 
The absorption imaging system consists of a CCD cam- 
era, a narrowband commercial diode laser at 1083 nm 
(TOPTICA, model DLIOO) and an AOM. The laser and 
AOM are used to generate low-intensity (/ < 0.05/sat) 
probe Hght pulses with a duration of 100 /xs, while the 
CCD camera (mounted behind a magnifier lens) records 
absorption images of the sample with a magnification 
of 0.17. The images are used to determine the (Gaus- 
sian) density distribution of the atoms in the trap |50|. 



n{x,y,z,t) = no{t) cxp 



z 

2^ 



(35) 



with no(t) the (time-dependent) density in the center 
{x = y = z — 0) of the sample, and Up and CTz the rms 
radii of the distribution; the number of trapped atoms 
N^JJJn{r)d\^i2nf/^ala, no. 

The MCP detectors allow for the independent monitor- 
ing of ions and He* atoms that escape the trap. With 
an exposed negative high voltage on its front plate, one 
MCP detector attracts all ions produced in the trapped 
sample. The other MCP is mounted behind a grounded 
grid and detects only the He* atoms that exit the trap in 
its direction. This shielded MCP detector is used to deter- 
mine the temperature T of the trapped samples through 
time-of-flight measurements. The unshielded MCP de- 
tector measures the instantaneous ionization rate in the 
trapped sample and is used to determine trap loss and 
ionization rates in the sample. 



A. Trap loss and ionization in magneto-optical trap 

For magneto-optically trapped He* samples, the time- 
evolution of the number of trapped atoms N can be de- 
scribed by the phenomenological equation 0, 



dN{t) 
dt 



L-aN{t)- (3 / / /n2(r,i)dV, 



(36) 



Writing down this differential equation, we assume that 
N{t) is controlled by three simultaneously occurring pro- 
cesses, corresponding to the three terms on the right- 
hand side of the equation. The first term is a constant 
loading rate L, representing the capture of atoms from 
the decelerated atomic beam into the mot. The sec- 
ond and third term are the linear and quadratic trap 
loss rate, respectively. The loss processes are defined in 
terms of the local atomic density of the sample n{r,t) 
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(with L = 0) by 



dn(r, t) 



-an{r,t) — (3n'^{r,t). 



(37) 



The nomenclature of the loss rate terms refers to their 
proportionality to density n{r,t) and density squared 
in?[r,t). Analogously, the proportionality constants, a 
and /3, are referred to as the linear and quadratic loss 
rate coefficient, respectively. 

For He* samples in a 1083 nm mot, only coUisional loss 
mechanisms give rise to significant trap loss. Quadratic 
trap loss is determined by collisions between trapped He* 
atoms, while linear trap loss results from collisions with 
particles traversing the trapping volume, such as back- 
ground gas particles and helium atoms from the atomic 
beam. Loss rates in ''He* samples can be estimated us- 
ing cross section data reported in literature [23. Is^ . Ta- 
bic IIVI presents an overview of trap loss mechanisms in 
■*He* samples; cross section data are combined with the 
atomic density in the sample, the background gas density 
or the atomic beam flux to determine the estimates. 

The trap loss mechanisms can be subdivided into 
ionizing mechanisms and mechanisms where atoms are 
lost without the formation of ions. Linear trap loss is 
dominated by a non-ionizing loss mechanism. As the 
beam of metastable atoms is not separated from the 
beam of ground state atoms (contrary to other work 
[lol IT3L HI . I2H'). collisions of ground state ''He atoms 
from the atomic beam with trapped atoms (collision en- 
ergy K, 4.9 meV) give rise to a trap loss rate of about 
2 s~'. Ionizing mechanisms hardly contribute: collisions 
with slowed {E^ k, 0.0064 meV) and non-slowed ^He* 
atoms {E-c w 6.5 meV) from the atomic beam, and colli- 
sions with background molecules (presumably H2O, H2, 
N2 and O2) result in a loss rate of about 1 x 10^^ s~'. 

Quadratic trap loss is dominated by ionization mech- 
anisms. In the presence of trapping light, light-assisted 
ionizing collisions between trapped *He* atoms give rise 
to an ionization rate of 4 s~'. In the absence of trapping 
light, ionizing collisions yield a rate of 0.1 s~^. As two 
atoms are lost for every ion that is formed, the corre- 
sponding trap loss rates are 8 s~^ and 0.2 s^', respec- 
tively. The escape of fast "^Hc* atoms from the trap, 
through fine-structure-changing collisions or radiative es- 
cape, constitutes a non-ionizing quadratic loss mecha- 
nism in the presence of trapping light, that can be ne- 
glected IsHssf . 

The ionization rate associated with the linear and 
quadratic ionization mechanisms can be written 



mon{t) ^ eh/3 



r?{r,t)^r, (38) 



where Ca and are the weights of ionization mecha- 
nisms in linear and quadratic trap loss, respectively. The 
quadratic ionization rate is half of the ionizing quadratic 
trap loss rate, as a single ion is formed for every pair of 
colliding He* atoms that is lost from the trap. As the 



linear ionization rate is small compared to the quadratic 
ionization rate, both in the presence and absence of trap- 
ping light, and quadratic trap loss is almost completely 
determined by ionization mechanisms, we can set, to 
good approximation, = and = 1, so that 



dMon(t) 

di 



(39) 



For samples of ^He* atoms, trap loss and ionization 
are determined by the same loss mechanisms. As collision 
studies are rare for "^He* atoms 53], we cannot derive 
accurate estimates of trap loss and ionization rates in this 
case. However, cross sections are not expected to show 
a large isotopic dependence (less than a factor of two), 
so we conclude that Eq. H39I) also applies to ^He* sam- 
ples. It should be noted in this respect, that hyperfine- 
changing collisions are forbidden by energy conservation. 
Trapped '^He* atoms are in the lower / = | hyperfine 
level, so that hyperfine-changing collisions require an en- 
ergy larger than the hyperfine splitting i?hf = 28 /xeV, 
which corresponds to a temperature iShf/f^B = 0.2 K. 



B. Ionization rates for light-assisted collisions 

To determine the ionization rate coefficient for light- 
assisted collisions, £6/3/2 « /3/2, we have performed a 
trap loss experiment, where the loading of atoms to the 
MOT is abruptly stopped and the decaying ionization rate 
in the sample is monitored with the unshielded MCP de- 
tector. For Gaussian samples, described by Eq. (|35|) . the 
ionization rate can be written as 



djVion(t) 

At 



^ V 



JL 
4^/2 



(40) 



where effective volume V = {2iTY^^apaz, such that cen- 
tral density no(i) — N{t) /V. The current signal from the 
MCP detector is proportional to the ionization rate and 
is converted to a voltage signal that is given by 



(41) 



where e is the electron charge and Rcs is an effective 
resistance. 

Substitution of Eq. H35|l into Eq. (|36|l shows that the 
central density satisfies differential equation 



dno{t) _ L _ /3 2 



dt V 



ano{t) 



2^/2 



(42) 



If the loading of atoms to the mot is abruptly stopped, 
the density decays as 



no{t) = 



no(0) 



(3no{0) \ /3no(0) 
1 H ^ — I exp(ai) 



2V20 



(43) 



2V20 
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TABLE IV: Trap loss rate dN/dt and ionization rate dA'ion/dt associated with collisional loss mechanisms in magneto-optically 
trapped samples of ''He* atoms. Rates are calculated for collisions of trapped ''He* atoms with ground state (1 ^S) 'He atoms 
and metastable (2 "^S) ''He* atoms from the atomic beam, for collisions with ground state ''He atoms and several molecules from 
the thermal background gas, and for collisions between trapped *He* atoms from the sample, both in the presence and absence 
of trapping light (with a wavelength of 1083 nm). For each collisional mechanism, collision energy Er is given, as well as total 
cross section a^*"*) or ionization cross section ct*^'""'. Furthermore, beam flux F, background gas density h, or central density 
no of the sample are given. For collisions with *He* atoms from the beam, we distinguish between slowed {Er = 0.0064 meV) 
and non-slowed {Er = 6.5 meV) atoms [soll . The cross sections are taken from various references; in references where an 
ionization rate coefflcient K is given the corresponding cross section is calculated as ct^'™-* = K/v^, where Vr = {fik-sT /■k^)^^'^ 
is the mean relative velocity for collisions of particles with mass mi and m2 in a gas with temperature T, with reduced mass 
^ = mim2/(mi + m2). 



Collision 


Originating 


Er 


^(tot) 




F 


n 


no 


diVion/dt 


dN/dt 


partner 


from 


(meV) 


(lO-^'^cm^) 


(10-'«cm2) 


(cm-^s-') 


(cm-3) 


(cm-^*) 






^He 


beam 


4.9 


200" 




-10" 









-2" 


*He* 


beam 


6.5 




ISl'^ 


4 X 10" 






7 X 10~^ ' 


1 X 10"^ 


^He* 


beam 


0.0064 




1160' 


8 X 10*^ 






9 X lO-'" 


2 X 10"^ 


*He 


background 


16.5 


140-^ 






1 X lO'' 







0.02» 


Ha 


background 


22 




0.1'' 




< 1.2 X lO*' 




< 2 X 10"'" 


< 4 X 10"'' 


H2O 


background 


6.0 




131' 




< 1.2 X 10^ 




< 9 X 10"''" 


< 2 X 10"^ 


O2 


background 


4.7 








< 1.2 X 10*^ 




< 4 X 10"^ s 


< 8 X 10"^ "* 


N2 


background 


4.1 




2*= 




< 1.2 X lO*' 




< 1 X 10"^ " 


< 2 X 10"^ 


^He* 


sample 


0.00001 




4 X 10^ ' 






4 X 10^ 


0.1" 


0.2'* 


''He* + light 


sample 


0.00001 




16 X 10'' ' 






4 X 10^ 


4'" 





""Vrincoanu et al. l56l . 

''For collisions with atoms from the beam, rates are aF. 

"Yor collision energies 0.1 meV < Er < 100 meV, erf'"") satisfies 
the (semi-classical) energy dependence cr('™) = ct i {E \/Er)" , with 
El = 1 meV, 0-1 = 318 X 10"^^ cm^ and a = 0.3 Is^. 

■^As two atoms are lost from the trap for each ion that is formed, 
dN/dt = 2{dN;^Jdt). 

"^Venturi et al. vm . 

•^Mastwijk ^ and Rothe et al. IHl . 

^For collisions with background gas particles, rates are aVrfi, 
where = {2Er / fi)^^'^ is the mean relative velocity for collisions 
of particles with mass m\ and m2, with ^ = mi 7712 /(mi -I- ra2) 
their reduced mass. 

'•Martin et al. Isl . 

'Mastwijk l57l : the large cross section results from a large attrac- 
tive force between ''He* and H2O, a consequence of the permanent 
dipole mome nt o f the H2O molecule. 

jParr et al. IgcII . 

*Yamazaki et a l. I61l . 

'ToI et al. lid : the cross section for light-assisted collisions is 
exceptionally large due to the optical excitation of the colliding 
atom pair to long-range dipole-dipole potentials 1^. 
"For samples with a Gaussian density distribution, as 
given by Eq. I35i . quadratic ionization rate dAfjon/dt = 
((T('°'i)t;r/Af) /// n2(r) d^r = crt'"") Vi.no/2\/2. 



Substitution of Eq. (|43|) into Eq. (|41|) gives an expres- 
sion for the decaying ionization signal. Loss rates a and 
/3no(0) determine the exact shape of the decay and are 
derived from the measured traces by means of curve fit- 
ting. 

We have performed trap loss measurements on '^He* 
and ^He* samples. The loading of atoms into the mot 
is stopped by blocking the Zeeman slowing light with 
the AOM used for frequency detuning the slowing light 
from the atomic transition. Decay traces are averaged 



four times using a digital oscilloscope. It has been ver- 
ified that the variations in the central density are small 
enough that an averaged decay curve allows an accurate 
determination of loss rates a and /3no(0). An example 
of an averaged decay trace and the corresponding fit are 
displayed in Figure[51 The central density in the samples 
no(0) is derived from absorption imaging. The resulting 
rate coefficients a and [3 are presented in Table Ivl 

The linear loss rate coefficients are close to 1 s~^ 
for both isotopes. Assuming that the loss rate stems 
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FIG. 5: Ionization signal for a trap loss measurement on a 
^He* sample. At t = 0, the input of atoms to the MOT is 
stopped abruptly. The rapid decay of the ionization signal 
(black dots) is non-exponential, as can be seen in the inset. 
Fitting the signal to our model (grey line) yields trap loss 
rates /3no(0) = 20 s"^ and a = 0.7 s~\ 



TABLE V: Characteristic parameters of the magneto- 
optically trapped He* samples and extracted loss rate co- 
efficients. Experimental errors correspond to one standard 
deviation. 



^Hc* 



*He^ 



T (mK) 
N{-) 
no (cm^^) 
a (s-i) 
/3 (cmVs) 
K (cmVs) 

7^("npol) (cj^Vs) 



2.0(3) 
2.6(9) X 10** 
3.0(5) X 10^ 

0.8(2) 
5.5(8) X 10"'' 
1.8(3) X 10"^° 
1.9(3) X 10"^° 



1.9(1) 
3.7(5) X 10® 
4.4(4) X 10** 

0.6(3) 
3.3(7) X 10"^ 
8(2) X 10"" 
10(2) X 10~" 



from collisions with ground state atoms from the atomic 

beam, we can determine the intensity of the beam of 

ground state atoms. Using the total cross section a = 

200 X 10~^^ cm^ js^ for collisions between trapped "^He* 

atoms and ground state *He atoms from the atomic 

beam, the distance between source and trapped sample 

of 370 cm, and loss rate coefficient a = 0.6 s~-^, we de- 
ls 



rive an intensity of 4 x 10 s sr 



As the intensity of 

^He* atoms is 4 x lO" s'^ sr"! [S^, the fraction of ^He* 
atoms in the beam is 10~^. 

The quadratic loss rate coefficient for '^He* is almost 
twice as large as the loss rate coefficient for ^He*. It has 
been pointed out that this isotopic difference stems from 
a difference in the relative number of ionization chan- 
nels, opened up by the lowering of centrifugal barriers by 
the long-range, resonant dipole-dipole interaction [ll| . 
Quantum statistical symmetry requirements play a role 
in these collisions, but the effects are washed out as the 
number of participating partial waves is much larger than 
one. The loss rate coefficient for ^He* is in good agree- 



ment with other work p^. 



C. Ionization in the absence of trapping light 

Collisions where ionization is not preceded by absorp- 
tion of trapping light contribute little to the ionization 
rate in magneto-optically trapped He* samples. How- 
ever, the corresponding ionization rate coefficients can 
be determined from a comparative measurement of the 
ionization rate in the presence and absence of trapping 
light 0. In this measurement, the trapping and slow- 
ing light are blocked for a short time interval of 100 /iS 
with the AOMs used for the detuning of the laser frequen- 
cies from the atomic transition. With the trapping light 
present, the observed ionization signal 



Von{t) = eR, 



eff 



4^2 



(44) 



is dominated by light-assisted collisions and is relatively 
large. If the trapping light is absent, collisions of '^He* 
or ^He* atoms occur without optical excitation; the cor- 
responding ionization rate is much smaller and can be 
written as 



K 
2V2 



4(0) 



(45) 



with K the ionization rate coefficient in the absence of 
trapping light. Combing Eqs. (|44|l and (|45|) . the rate 
coefficient can be written as 



2 ipo 



(46) 



Clearly, K can be derived from a measurement of the ion- 
ization rate coefficient for light-assisted collisions, /3/2, 
and signals (fon, ^oS and (phgr- 

In the experiment, the trapping and slowing light are 
blocked every 200 ms, and the ionization signal is aver- 
aged 256 times with a digital oscilloscope. Examples of 
averaged ionization signals are displayed in FigureEl The 
measurement is repeated with the atomic beam blocked 
to obtain the background ionization signal fbgr, that in- 
cludes an offset originating from the MCP signal amplifier. 
To obtain accurate values of {ifon — 'fihgT) and (Voff— ^bgr), 
the average of the signals is determined over 40 fis time 
intervals, as indicated in Figure 

Although the atoms are not confined if the trapping 
light is absent, the expansion of the trapped He* sample 
during 100 /is is insignificant and can be neglected. It 
is checked that the trapped sample remains unaffected if 
the light is blocked repeatedly, running the experiment 
with a duty cycle of 200 ms. The resulting ionization rate 
coefficients are presented in Table The coefficients 
are close to but do not agree within the error bars with 
other experimental results ITollTTll2ll | . that suffer from 
mutual inconsistency themselves, as has been pointed out 

M. 
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FIG. 6: Ionization rates in the presence and absence of 
trapping (and slowing) light. The signals are averaged over 
40 /is time intervals indicated by the shaded areas to obtain 

{(fon - ^bgi) and (ifoS ~ V'bgr). 



D. Comparison between experiment and theory 

To confront the theoretical model of Sec. ^1 with the 
measurements of Sees. IIII Bl and IIII CI we compare the 
theoretical expression for K, Eq. H32|l . with the T — 2 mK 
experimental values of Table Ivl As optical pumping pro- 
cesses cause the magnetic substate distribution to deviate 
from a uniform distribution (with the stretched substates 
slightly overpopulated), we need to determine the mag- 
netic substate distribution Pm(r) in the trapped samples 
to calculate a value for K that can be compared to ex- 
periment. 

Distribution P„i{r) is obtained as the steady-state so- 
lution of a rate equation model, that describes the opti- 
cal pumping in the MOT |2lj | . Starting from the intensity 
and detuning of the trapping light, and an expression 
for the quadrupole magnetic field, rate equations are for- 
mulated and subsequently solved to obtain the steady 
state substate population in the sample. In these cal- 
culations we take into account that the intensities of the 
trapping beams are not balanced and that (consequently) 
the trapped sample is not exactly centered on the zero 
point of the magnetic field. 

At a temperature T — 2 mK, the resulting theoretical 
ionization rate coefficients. 



KsHc = 1.7 X 10-1° cmVs, 
KiHc = 8.0 X 10"" cmVs, 



(47) 
(48) 



agree very well with the experimental values of Table Ivl 
At this low temperature only s and p-waves have to be 
taken into account as higher order waves contribute less 
than a few percent. The samples are partially polarized, 
so that quadratic ionization rates are smaller than they 
would be in unpolarized samples. Because of the different 
substate structure (s = 1 compared to / = 3/2) and 
different quantum statistical symmetry, the corrections 



for the two isotopes differ, i^4Ho/-f^4Hc 
KsHc/ K^^f°^^ = 0.93. For direct comparison with the 
theoretical results of TableEIlX^""?"') for ^He* and ^He* 
have been included in Table IVl as well. 



IV. DISCUSSION AND CONCLUSIONS 

The theoretical model presented in Sec. ITTI shows good 
agreement with other theoretical work and with the ex- 
perimental results reported in Sec. IIIII The agreement 
between the various results has to be considered partly 
coincidental in view of the theoretical uncertainties in 
the potentials, as already mentioned in Sec. IIIFI Any- 
way, this shows that cold ionizing colhsions of He* atoms 
can be understood as single-channel processes that are 
determined by Wigner's spin-conservation rule, quantum 
threshold behavior and the symmetrization postulate. 
Using the model, the difference between the ionization 
rate coefficients for "^He* and *He* can be interpreted as 
an effect of the different quantum statistical symmetry 
of the two isotopes and the presence of a nuclear spin 
in the case of •^He*; these differences do not depend on 
uncertainties related to the potential as the same po- 
tentials (with mass scaling) were used in the ^He* and 
■^He* case. As the model is relatively simple, it is com- 
plementary to the more complete (and precise) close- 
coupling theor y th at has been developed for ^He* col- 
lisions [idL Ha. Il7l|. Another single-channel model, that 
was published recently was shown to be based on 
erroneous assumptions. 

The experimental values of Sec. Illll do not agree 
with other (mutually inconsistent) experimental results 
[^ ITollllll2H . The discrepancy between the experimental 
values is difficult to interpret, as the experimental deter- 
mination of the ionization rate coefficient is the result 
of an extensive analysis, including the determination of 
the density distribution and temperature of the trapped 
samples, as well as the distribution over the magnetic 
substates. The ionization rate coefficient is particularly 
sensitive to the density in the sample. It must be checked 
carefully if frequency drifts of the probe laser light or 
stray magnetic fields are small enough to avoid incom- 
plete absorption, which would result in underestimation 
of the density and overestimation of the ionization rate 
coefficients. If the number of trapped atoms is obtained 
from fluorescence imaging |ll| . an accurate calibration 
of the CCD chip is crucial. Furthermore, collisions of 
trapped atoms with background atoms or atoms from the 
atomic beam must be considered. If the trapped atom 
number is small (< 10^), quadratic ionization becomes 
small and other (linear) ionization mechanisms possibly 
play a part, hampering accurate measurements. Finally, 
it should be noted that in most experiments d, 0, ^ 
the magnetic substate distribution has not been taken 
into account. 

It would be interesting to extend the work presented 
here to trapped samples containing both isotopes and 
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study heteronuclear ionizing collisions. In the case of 
collisions between distinguishable particles, quantum sta- 
tistical symmetry requirements should be absent, which 
could be confirmed from an investigation of ionizing col- 
lisions. This work is in progress in our laboratory. An- 
other interesting extension of the work presented here, is 
the study of the ionization rates for samples with a pre- 
pared substate population. It raight be possible to study 
depolarization due to collisions [13 ■ Finally, ionizing col- 
lisions can also be investigated in the quantum degener- 
ate regime |62- Ionization rates could be used to study 
quantum statistical properties of a quantum degenerate 
mixture with high spatial and temporal resolution. It is 



conceivable that phase separation in a quantum degen- 
erate mixture will be observable through the ionization 
rate in the sample. 
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